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Abstract 

Recently it was shown that there is a unique 'L^ symmetry for the MSSM which 
ahows the Yukawa coupHngs and dimension five neutrino mass operator, forbids 
the n term and commutes with SO(IO). This 'L^ symmetry contains matter parity 
as a subgroup and forbids dimension four and five proton decay operators. We 
show how to construct string vacua with discrete R symmetries in general and 
this symmetry in particular, and present an explicit example which exhibits the 
exact MSSM spectrum, the 'L^ symmetry as well as other desired features such as 
gauge-top unification. We introduce the Hilbert basis method for determining all 
D-flat configurations and efficient algorithms for identifying field configurations 
with a desired residual symmetry. These methods are used in an explicit example, 
in which we describe in detail how to construct a supersymmetric vacuum configu- 
ration with the phenomenologically attractive 'L^ symmetry. At the perturbative 
level, this is a supersymmetric Minkowski vacuum in which almost all singlet fields 
(moduli) are fixed. 



1 Introduction 



There are many independent observations hinting at the relevance of a high scale for par- 
ticle physics. The smallness of neutrino masses has a simple explanation in terms of the 
see-saw mechanism pT] and relies on the existence of heavy singlet neutrinos. Stabilizing 
the electroweak scale against the see-saw scale seems to require supersymmetry; remark- 
ably, the simplest supersymmetric extension of the standard model, the MSSM, realizes 
the compelling scenario of gauge unification [2] at a scale Mqut — 2 ■ 10^^ GeV, which is 
suspiciously close to the see-saw scale. Both scales are not too far from Mp ~ 2-10^'^ GeV, 
which is set by Newton's constant. 

The question of how to incorporate all scales in a coherent scheme has been addressed 
for more than 30 years. From a bottom-up perspective one is led to the scheme of grand 
unified theories (GUTs). Although this scheme exhibits various very appealing features, 
there are three major obstacles. First, there is the so-called doublet-triplet splitting, 
and related to it, the MSSM /i problem. Second, even if this problem is solved, unified 
models typically are in conflict with dimension five proton decay ^H] operators, fit 
is well known that dimension four proton decay can be forbidden by matter paritylj) 



Third, in four- dimensional models of grand unification there is no relation between the 
GUT and Planck scales, Mqut and Mp. 

String theory is believed to provide us with such a relation. However, if string 
theory is to describe the real world, it should also provide us with solutions to the first 
and second problems. In fact, as known for a long time, the doublet-triplet splitting 
problem has a simple solution in theories with extra dimensions in which the GUT 
symmetry is broken in the process of compactification [51E], which also avoids the most 
stringent problems with dimension five proton decay [7]. However, in concrete string 
compactifications (see P4T0] for early attempts and [II1|12] for a different approach) 
very often the problem is reintroduced; this applies also to the models discussed more 
recently [T3HT6]. On the other hand, one cannot rule out these constructions as their 
vacua are not completely understood. That is, the analysis of potentially realistic string 
models is a non-trivial task since a given model exhibits a plethora of vacua with very 
different features. The role of discrete symmetries in identifying and analyzing such 
vacua has been stressed recently [TTj. One of these symmetries is matter parity, which has 
been successfully embedded in string theory [16]. This study is devoted to a discussion 
of the role of further discrete symmetries in such models and their phenomenological 
implications. Specifically, we will focus on a proposed "L^ symmetry [TS] which has 
recently been shown in to be the unique anomaly-free possibility with the following 
properties: 

1. it forbids the /i term at the perturbative level; 

2. it allows the MSSM Yukawa couplings and the effective neutrino mass operator; 

^Matter parity is sometimes also known as "i? parity". We choose not to use this terminology as 
matter parity is non-i?, i.e. it commutes with supersymmetry, and one point of this paper is to discuss 
a true discrete i? symmetry. 
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3. it commutes with SO(IO) in the matter sector. 

This symmmetry has the appeahng feature that it forbids automatically dimension four 
and five proton decay operators. We will discuss how to identify string vacua exhibiting 
this symmetry and present a globally consistent string-derived model with the exact 
MSSM spectrum realizing this symmetry. 

We start in section [2] with a short description of the general picture. In section 12] we 
present an explicit string- derived model in which an anomalous symmetry explains a 
suppressed vacuum expectation value of the superpotential, provides us with a solution of 
the n problem and suppresses dimension five proton decay operators. Section H] contains 
our conclusions. In various appendices we collect details of our calculations. 

2 General picture 

String theory compactifications provide us with a plethora of vacuum configurations, 
each of which comes with symmetries and, as a consequence, with extra massless degrees 
of freedom whose mass terms are prohibited by these symmetries. Simple examples for 
such compactifications include heterotic orbifolds [201121], where the rank of the gauge 
group after compactification equals that of Eg x Eg, i.e. 16. A few hundreds of orbifold 
models are known in which Eg x Eg gets broken to the standard model gauge symmetry 
= SU(3)c X SU(2)l x U(l)y (with hypercharge in GUT normalization) times U(l)" 
times a hidden sector group and the chiral spectra of the MSSM [T51I22]- They also 
exhibit exotics which are vector-like with respect to Gsm and which can be decoupled 
when the extra gauge symmetries are broken. Each of these models contains many 
vacua, i.e. solutions of the supersymmetry conditions Vp = Vd = 0. Typically these 
vacua exhibit flat directions before supersymmetry breaking. 

At the orbifold point, where the vacuum expectation values (VEVs) of all fields are 
zero, we have discrete R as well as continuous and discrete non-i? symmetries. Typically 
one of the U(l) symmetries appears anomalous, which is conventionally denoted by 
U(l)anom- Also some of the discrete symmetries may appear anomalous [23l|2l]. After 
assigning VEVs to certain fields, some of the symmetries are spontaneously broken and 
others remain. We shall be mainly interested in remnant discrete symmetries, which 
can be of R or non-i? type and be either anomalous or non-anomalous. We will discuss 
examples of all kinds in section [31 

Clearly, one cannot assign VEVs to the fields at will. Rather, one has to identify 
field configurations which correspond to local minima of the (effective) scalar potential. 
Let us briefly describe the first steps towards identifying such vacuum configurations. 
Consider a configuration in which several fields attain VEVs. We focus on "maximal 
vacua" (as in [I7|), i.e. we assume that all fields which are neutral under the remnant 
gauge and discrete symmetries, called (1 < < A^) in what follows, attain VEVs 
(if these are consistent with D-fiatness). All fields without expectation value, denoted 
by (1 < J < M), therefore transform non-trivially under some of the remnant 
symmetries. 
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2.1 Discrete non-i? symmetries 

The case of vacua with non-R discrete symmetries has been discussed in detail in [T^[T7] . 
In this case, the superpotential has the form 

W = . . . + (terms at least quadratic in the ip^^^) . (2.1) 

Therefore, the F-term equations for the ^'^^^ fields trivially vanish and we are left with 
F-term equations for the 0*^*^ fields, which generically have solutions. Hence, if all (p^^^ 
enter gauge invariant monomials composed of (f)^'^^ fields only, we will find supersymmetric 
vacua, i.e. solutions to the F- and D-term equations. 

Because of the above arguments it is sufficient to look at the system of 0*^*^ fields only, 
which has been studied in the literature. Consider the case of a generic superpotential 
W. It is known that the solutions to the D- and F-term equations intersect generically 
in a point [25]. That is, there are point-like field configurations which satisfy 

Da = = at 0» = , (2.2) 

where, as usual. 

Da = 5^(0«)*T,0«, (2.3a) 
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The term 'point-like' means that there are no massless deformations of the vacuum (12. 2p . 
The reason why these vacua are point-like is easily understood: generically the F-term 
equations constitute as many gauge invariant constraints as there are gauge invariant 
variables. However, this also means that, at least generically. 



^ . (2.4) 



If the fields attain VEVs (</)*-*^) of the order of the fundamental scale, one hence expects to 
have too large a VEV for W. One possible solution to the problem relies on approximate 
R symmetries [SSj, where one obtains a highly suppressed VEV of the superpotential. 
In what follows, we discuss an alternative: in settings with a residual R symmetry the 
above conclusion can be avoided as well. 

2.2 Discrete R symmetries 

Let us now discuss vacua with discrete R symmetries. To be specific, consider the order 
four symmetry Tj^, under which the superpotential W has charge 2, such that 

W ^ -W (2.5) 
under the Tjf generator (. Superspace coordinates transform as 

9a iOa (2.6) 
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such that the F-term Lagrangean 

= jd^eW + h.c. (2.7) 

is invariant. Chiral superfields will have R charges 0, 1, 2, 30 Both the fields of the type 
ipi and ip3 with R charges 1 and 3, respectively, can acquire mass as the ijjf and '03 terms 
have R charge 2 mod 4 and thus denote allowed superpotential terms. 

The system of fields ^q"* and ip2^ with R charges and 2, respectively, is more 
interesting. Consider first only one field 0o and one field ip2- The structure of the 
superpotential is 

^ = 02 ■ /(0o) + 0{^l) (2.8) 
with some function /. The F-term for (f)Q vanishes trivially as long as is unbroken, 
dW 

jr- = ^2- f'{<Po) = . (2.9) 

Note that due to the Tj^ symmetry the superpotential vanishes in the vacuum. Thus 
it is sufficient to look at the global supersymmetry F-terms. On the other hand, the 
F-term constraint (at 02 = 0) 

will in general fix 0o at some non-trivial zero (0o) of /. Indeed, there will be a super- 
symmetric mass term, which can be seen by expanding 0o around its VEV, i.e. inserting 
00 = (0o) + 50o into dlSD, 

^ = /'((0o)) 500 ^2 + 0(50^, 02')- (2.11) 

The supersymmetric mass /'((0o)) is generically different from 0. 

Repeating this analysis for 0q ■* and M ip^"^ fields reveals that the F-terms of the 
'02"''' lead to M, in general independent, constraints on the 0o'' VEVs. For N = M we 
therefore expect point-like vacua with all directions fixed in a supersymmetric way. 

To summarize, systems with a residual R symmetry ensure, unlike in the case without 
residual symmetries, that (W) = 0. However, in systems which exhibit a linearly realized 
Tjf somewhere in field space it may not be possible to find a supersymmetric vacuum that 
preserves Tj^. In the case of a generic superpotential this happens if there are more, i.e. 
M > N, fields with R charge 2 than with 0. On the other hand, if there are more fields 
with R charge than with 2, i.e. for M < N, one expects to have a Minkowski vacuum 
with N — M fiat directions. For N = M one can have supersymmetric Minkowski vacua 
with all directions fixed in a supersymmetric way. 

An important comment in this context concerns the moduli-dependence of couplings. 
As we have seen, in the case of discrete i?-symmetries one might obtain more constraint 
(i.e. F-term) equations than i?-even 'matter' fields. Specifically, in string vacua one 
should, however, carefully take into account all i?-even fields, also the Kahler and com- 
plex structure moduli, Tj and Uj, on whose values the coupling strengths depend. 



special role is played by the dilaton S, whose imaginary part a = lmS\e=o shifts under 
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3 An explicit string-derived model 



In order to render our discussion more specific, we base our analysis on a concrete model. 
We consider a Z2 x Z2 orbifold compactification with an additional freely acting of 
the Eg X Eg heterotic string. Details of the model including shift vectors and Wilson 
lines can be found in appendix [El 

In [27] a vacuum configuration of a very similar Z2 x Z2 model with matter parity 
and other desirable features was presented. However, the vacuum configuration discussed 
there has the unpleasant property that, at least generically, all Higgs fields attain large 
masses. In what follows we discuss how this can be avoided by identifying vacuum 
configurations with enhanced symmetries. In [19] another vacuum with the symmetry 
discussed in the introduction was found by using the methods presented in this paper. In 
both models the GUT symmetry is broken non-locally. This may be advantageous from 
the point of view of precision gauge unification [2S]. It also avoids fractionally charged 
exotics, which appear in many other compactifications (cf. the discussion in [29]). 

Labeling of states. We start our discussion with a comment on our notation. In 
a first step, we label the fields according to their Gsm x [SU(3) x SU(2) x SU(2)]hid 
quantum numbers. In particular, we denote the standard model representations with 
lepton/Higgs and d-quark quantum numbers as 



In the next step we identify such that the Li/Li decompose in lepton doublets 
with odd charges and Higgs candidates hd/hu with even charges etc. The details 
of labeling states are given in appendix |El 

Searching for How can one obtain vacua with in practice? We found the 

following strategy most efficient: 

1. In a first step we switch on a random sample of SM singlets in such a way that all 
unwanted gauge factors are spontaneously broken. 

2. With these VEVs at hand, the original gauge and discrete symmetries at the 
orbifold point get broken to a discrete subgroup, which can be determined unam- 
biguously with the methods described in [30]. Details of the automatization of 
these methods are explained in [3T] . 

3. We only keep configurations in which there is a residual symmetry with pre- 
cisely three generations of matter having i?-charge 1; details of how to identify 
such configurations are given in appendix [Bi 



(l,2)-l/2, 

(1,2)1/2, 

(3,1)^1/3, 

(3,1)1/3. 



(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 
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4. From these configurations we select those exhibiting the following properties: 

• F- and D-fiat; 

• all exotics decouple; 

• one pair of massless Higgs, i.e. /x term forbidden to all orders (at the pertur- 
bative level); 

• Yukawa couplings have full rank. 

One of the main achievements of this study is a considerable simplification in the verifica- 
tion of the four items listed in step 4. In order to check D-fiatness of a given configuration 
we use the Hilbert basis method, which is described in detail in appendix O The other 
three properties can be verified by inspecting the remnant discrete symmetries only. In 
earlier studies [T314T6] we had to explicitly identify couplings that are consistent with 
the string selection rules in order to show that all exotics decouple and the Yukawa 
couplings have full rank. In our new approach the remnant symmetries will tell us im- 
mediately whether an entry of a mass or Yukawa matrix will or will not appear. We 
have cross-checked this method extensively by explicitly computing the couplings be- 
tween the charged and the VEV fields, and were always able to find a coupling which 
fills in an entry of a matrix, albeit sometimes at very high orders. Note, we assume that 
all couplings allowed by string selection rules appear in the superpotential. 

VEV configuration. Following the above steps, we obtained a promising configura- 
tion in which the fields 

0*'*'' = {01,02,03,04,05,06,07,08,09,010,011,012,013,014, 

Xl, X2, X3, X4, X5, Xi, X3, X4, X5, 1/3, 1/4, I/5, T/e} (3.2) 

attain VEVs. The full quantum numbers of these fields are given in table IE. 21 in ap- 
pendix [El In order to ensure D-flatness with respect to the hidden sector gauge factors, 
in a given basis not all components of the Xi/xi and attain VEVs. Details are given 
in equations (1D.3P and (ID.Sp in appendix iDl 

Remnant discrete symmetries. By giving VEVs to the fields in (13. 2p . we arrive 
at a vacuum in which, apart from Gsm and a 'hidden' SU(2), all gauge factors are 
spontaneously broken. The vacuum exhibits a symmetry, whereby the superpotential 
W has charge 2. 

The charges of the matter fields are shown in table [3TT1 The detailed origin of the 
Ijf symmetry is discussed later. Given these charges, we confirm by a straightforward 
field-theoretic calculation (cf. [2ll|32]) that appears indeed anomalous with universal 
SU(2)L-SU(2)L-Zf and SU(3)c-SU(3)c-Zf anomalies (see [IH] and appendix O]). 
The statement that appears anomalous means, as we shall discuss in detail below, 
that the anomalies are cancelled by a Green-Schwarz (OS) mechanism. On the other 
hand, the has a, by the traditional criteria, non-anomalous subgroup which is 
equivalent to matter parity [19] . 
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(a) Quarks and leptons. 



Qi di £{ Ci 

1 1 1 1 1 



(b) Higgs and exotics. 



hi h2 hs /i4 ha hi h2 h^ h4 h^ h^ 6i 62 S3 61 62 S3 

020200020022022200 



Table 3.1: charges of the (a) matter fields and (b) Higgs and exotics. The index i in 
(a) takes values i = 1, 2, 3. 

Z)-flatness. As already discussed, we cannot switch on the (p^'^^ fields at will; rather we 
have to show that there are vacuum configurations in which all these fields acquire VEVs. 
This requires to verify that the D- and F-term potentials vanish. With the Hilbert basis 
method (see appendix [C]) we could identify a complete set of D-fiat directions composed 
of 0^*-* fields. We compute the dimension of the Z)-fiat moduli space using Singular [53] 
and the STRINGVACUA [31] package; the result is that there are 18 Z)-fiat directions; 
the details of the computation are collected in appendix |Dl 

i^-term constraints. Next we consider the F-term constraints. As discussed in sec- 
tion 121 the F-term conditions come from the fields with i?-charge 2. We compute the 
number of independent conditions in appendix |Dl The result is that there are 23 in- 
dependent conditions on 18 + 6 = 24 D-fiat directions, where we included the Kahler 
and complex structure moduli. We therefore expect to find supersymmetric vacuum 
configurations in which all the (f)^^'' acquire VEVs. In this configuration, almost all sin- 
glet fields, including the geometric moduli are fixed in a supersymmetric way. It will 
be interesting to compare this result to similar results found recently in the context of 
smooth heterotic compactifications [35]. We expect a significantly different, i.e. health- 
ier, phenomenology than in the case in which a large number of singlets acquire mass 
only after supersymmetry breaking [36(137]. Notice that there are two possible caveats. 
First, the analysis performed strictly applies only to superpotentials which are, apart 
from all the symmetries we discuss, generic. Second, it might happen that there are 
supersymmetric vacua, but they occur at large VEVs of some of the fields, i.e. in regions 
of field space where we no longer control our construction. Both issues will be addressed 
elsewhere. 

Higgs vs. matter. The subgroup of the Tj^ symmetry allows us to discriminate 

between 

• 3 lepton doublets, ii = {L4^, Lq, Lj}, 
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• 3 d-tjpe quarks, di = {Di, D^, D4}, 
on the one hand, and 

• Higgs candidates, hi = {Li, L2, L3, L5, Lg, Lg} and hi = {Li, L2, 13, L4, L5, Lq}, 

• exotic triplets, 6i = {Di, D2, D^} and 6i = {D2, D^, D^} 
on the other hand. 



Decoupling of exotics. With the charges in table 13.11 we can readily analyze the 
structure of the mass matrices. We crosscheck these structures by explicitly computing 
the couplings allowed by the string selection rules (cf. |2Z])- Note there is a caveat: our 
results are based on the assumption that all couplings that are allowed by the selection 
rules will appear with a non- vanishing coefficient. A (f)^ in the matrices represents 
a known polynomial of order n in the fields which we have calculated using string 
selection rules. A zero entry in the matrices means that the corresponding coupling is 
not present in the perturbative superpotential. The hi — hj Higgs mass matrix is 



M, 



/O 

07 




03 














02 




h2 



\ 












Pl3 











(3.3) 



0' 08/ 



Here we omit coefficients, which depend on the three Kahler moduli Tj and complex 
structure moduli Ui. Clearly, this mass matrix has rank five, such that there is one 
massless Higgs pair 

hu = aihi + a2 hj, + a^h^ , (3.4a) 

hd = bihi + 62 hs + 63 /i5 + 64 he (3.4b) 
with tti and bj denoting coefficients. The 6 — 6 mass matrix is 

\ 












(3.5) 



Hence, the matrix has full rank and all exotics decouple. Note that the block structure 
of 5 is not a coincidence but a consequence of the fact that 62/63 and 62/63 form D4 
doublets (see below). Altogether we see that all exotics with Higgs quantum numbers, 
and all but one pair of exotic triplets, decouple at the linear level in the 0^*-* fields. This 
leads to the expectation that all but one pair of exotics get mass of the order of the 
GUT (or compactification) scale Mgut while one pair of triplets might be somewhat 
lighter. We also note that the presence of colored states somewhat below Mgut can give 
a better fit to MSSM gauge coupling unification (cf. [38] )• However, a crucial property 
of the 6- and 6 triplets is that, due to the Zf symmetry, they do not mediate dimension 
five proton decay. 
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Effective Yukawa couplings. The effective Yuk;awa couplings are defined by 



i=l,3,4 

The Yukawa coupling structures are 



yiX) 



i=l,3,5,( 



^6 \ 



/ 02 04 

y 06 06 
/ ^6 l6 26 \ 



y{3) 



a4 24 



k2 



06 

/ ;^6 



6 ^6 
f 



\ 



I \ 

06 

04 y 



(3.7a) 



(3.7b) 



(3.7c) 



Yd and Ye coincide at tree-level, i.e. they exhibit SU(5) GUT relations, originating from 
the non-local GUT breaking due to the freely acting Wilson line. There are additional 
contributions to Yu from couplings to /14 and to Y^/Yd from couplings to hi^^ which can 
be neglected if the VEVs of the 0''*-' fields are small. 

Because of the localization of the matter fields, we expect the renormalizable (1,3) 
and (3,1) entries in yJ^'* to be exponentially suppressed. 

Gauge-top unification. The (3, 3) entry of Y„ is related to the gauge coupling. More 
precisely, in an orbifold GUT limit in which the first orbifold plane is larger than 
the other dimensions there is an SU(6) bulk gauge symmetry, and the ingredients of the 
top Yukawa coupling hu (i.e. the fields /ii,3,4), uz and gs are bulk fields of this plane, 
i.e. hypermultiplets in the N = 2 supersymmetric description. As discussed in [39], this 
implies that the top Yukawa coupling yt and the unified gauge coupling g coincide at 
tree-level. Moreover, localization effects in the two larger dimensions [10] will lead to a 
slight reduction of the prediction of yt at the high scale such that realistic top masses 
can be obtained. 



flavor symmetry. The block structure of the Yukawa matrices is not a coincidence 
but a consequence of a D4 flavor symmetry [H], related to the vanishing Wilson line 
in the ei direction, Wi = (cf. e.g. |12])- The first two generations transform as a D4 
doublet, while the third generation is a 1^4 singlet. 



Neutrino masses. In our model we have 11 neutrinos, i.e. SM singlets whose charges 
are odd under Tjf meaning that they have odd charge, where is the matter 
parity subgroup of . Their mass matrix has rank 11 at the perturbative level. The 
neutrino Yukawa coupling is a 3 x 11 matrix and has full rank. Hence the neutrino 
see-saw mechanism with many neutrinos [l^ is at work. 
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Proton decay operators. The symmetry forbids all dimension four and five 
proton decay operators at the perturbative level [19]. In addition, the non-anomalous 
matter parity subgroup forbids all dimension four operators also non-perturbatively. 
The dimension five operators like qqq^ are generated non-pertubatively, as we will 
discuss below. 

Non-perturbative violation of Z^. Once we include the terms that are only for- 
bidden by the symmetry, we obtain further couplings. An example for such an 
additional term is the dimension five proton decay operator, 

where we suppressed coefficients. The bracket structure between the 0ii/0i2 and ^g/^io 
is a consequence of the non-Abelian D4 symmetry, where these fields transform as a 
doublet. The dot ' ■ ' indicates the standard scalar product. Note that there are in- 
variants with more than two D4, charged fields which cannot be written in terms of a 
scalar product. Further, S is the dilaton and the coefficient a = Svr^ in e~"'^ is such 
that e~° ^ has positive anomalous charge with respect to the normalized generator of the 
'anomalous' U(l). This generator is chosen such that it is the gauge embedding of the 
anomalous space group elemenlH (cf. equation flE.Sj) ). 

tanom = W3 + Eg X Eg lattice vectors . (3.9) 

The discrete Green-Schwarz mechanism is discussed in detail in |44] . 

Solution to the /x problem. The T^f anomaly has important consequences for the 
MSSM fi problem. The fi term is forbidden perturbatively by however, it appears at 
the non-perturbative level. Further, this model shares with the mini-landscape models 
the property that any allowed superpotential term can serve as an effective /i term (cf. the 
discussion in [SB] ) ■ This fact can be seen from higher-dimensional gauge invariance [IS] . 
Therefore, the (non-perturbative) /i term is of the order of the gravitino mass, 

;U ~ (^) ~ m3/2 (3.10) 

in Planck units. If some 'hidden' sector dynamics induces a non-trivial the /z 

problem is solved. 

In our model, we have only a 'toy' hidden sector with an unbroken SU(2) gauge group 
and one pair of massless doublets whose mass term is prohibited by Tjf. This sector has 
the structure discussed by Affieck, Dine and Seiberg (ADS) [16]. We find that the ADS 
superpotential is covariant. However, the hidden gauge group is probably too small 
for generating a realistic scale of supersymmetry breakdown. Yet there are alternative 
ways, such as the one described in [2S], for generating a hierarchically small 

■^See [24] for the discussion in a more general context. Note that we can always bring the anomalous 
space group element to the form (0'^w^,O) by redefining the model input appropriately. This amounts 
to a redefinition of the 'origin' of the orbifold. 
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Origin of In the orbifold CFT description tlie originates from the so-called 

if- momentum selection rules |17] (see also [HlllH]). These selection rules appear as 
discrete R symmetries in the effective field theory description of the model. We would 
like to stress that in large parts of the literature the order of these symmetries was given 
in an unfortunate way. This criticism also applies to the papers by some of the authors 
of this study. For instance, the orbifold plane was said to lead to a symmetry, but 
it turned out that there are states with half-integer charges. We find it more appropriate 
to call this symmetry Ijf, and to deal with integer charges only. In our model we have 
three Ijf symmetries at the orbifold point, stemming from the three Z2 orbifold planes. 

if-momentum corresponds to angular momentum in the compact space; therefore the 
discrete R symmetries can be thought of as discrete remnants of the Lorentz symmetry 
of internal dimensions. That is to say that the orbifold compactification breaks the 
Lorentz group of the tangent space to a discrete subgroup. In this study we content 
ourselves with the understanding that these symmetries appear in the CFT governing 
the correlators to which we match the couplings of our effective field theory. The precise 
geometric interpretation of this symmetry in field theory will be discussed elsewhere. 

The actual charges of [SU(3) x SU(2) x SU(2)]hid invariant expressions in this 
model are given by 

fe« = gx + -R2 + 2ra3, (3.11) 

where qx is the U(l) charge generated by 

tx = (4,0,10,-10,-10,-10,-10,-10) (-10,0,5,5,-5,15,-10,0) , (3.12) 

R2 denotes the R charge with respect to the second orbifold plane and is the lo- 
calization quantum number in the third torus. The relevant quantum numbers are 
given in table IE.2I The expression (13.111) for g^s is not unique, there are 17 linear 
combinations of U(l) charges and discrete quantum numbers which can be used to 
rewrite the formula without changing the T^f charges. Also the U(l) factors contained 
in [SU(3) X SU(2) x SU(2)]hid can be used to redefine tx- We refrain from spelling this 
out as we find it more convenient to work with invariant monomials (cf. the discussion 
in appendix [D]). It is straightforward to see that all monomials we switch on have R 
charge 0. 

4 Summary 

We have re-emphasized the important role of discrete symmetries in string model build- 
ing. As an application, we discussed an explicit string model which exhibits MSSM 
vacua with a symmetry, which has recently been shown to be the unique symmetry 
for the MSSM that forbids the /i term at the perturbative level, allows Yukawa couplings 
and neutrino masses, and commutes with SO(IO). This has a couple of appealing 
features. First, the n term and dangerous dimension five proton decay operators are 
forbidden at the perturbative level and appear only through (highly suppressed) non- 
perturbative effects. Second, at the perturbative level, the expectation value of the 
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superpotential is zero; a non-trivial expectation value is generated by non-perturbative 
effects. These two points imply that /i is of the order of the gravitino mass m3/2, which 
is set by the expectation value of the superpotential (in Planck units). 

The model is a x orbifold compactification of the Eg x Eg heterotic string. 
We discussed how to search for field configurations which preserve Tjf and how to find 
supersymmetric vacua within such configurations. The Hilbert basis method allowed 
us to construct a basis for all gauge invariant holomorphic monomials, and therefore to 
survey the possibilities of satisfying the D-teim constraints. As we have seen, in the case 
of residual R symmetries it may in principle happen that the F-term equations over- 
constrain the system. We have explicitly verified that this is not the case in our model, 
i.e. there are supersymmetric vacua with the exact MSSM spectrum and a residual l^f 
symmetry. Let us highlight the features of the model: 

• exact MSSM spectrum, i.e. no exotics; 

• almost all singlet fields/moduli are fixed in a supersymmetric way; 

• non-local GUT breaking, i.e. the model is consistent with MSSM precision gauge 
unification; 

• dimension four proton decay operators are completely absent as Tjf contains the 
usual matter parity as a subgroup; 

• dimension five proton decay operators only appear at the non-perturbative level 
and are completely harmless; 

• the gauge and top- Yukawa couplings coincide at tree level; 

• see-saw suppressed neutrino masses; 

• /i is related to the vacuum expectation value of the superpotential and therefore 
of the order of the gravitino mass; 

• there is an SU(5) GUT relation between the r and bottom masses. 

There are also two drawbacks: first, there are also SU(5) relations for the light gener- 
ations and second the hidden sector gauge group is only SU(2) and therefore probably 
too small for explaining an appropriate scale of dynamical supersymmetry breaking. 

Although we have obtained a quite promising string vacuum, the main focus of this 
paper was on developing new methods rather than working out the phenomenology of a 
model. We have discussed in detail how to determine the residual discrete symmetries 
of a given VEV configuration. As a consequence, we could immediately understand 
the features of such a configuration. For instance, in earlier studies [T3HT6] we had to 
explicitly identify couplings that are consistent with the string selection rules in order 
to show that all exotics decouple and the Yukawa couplings have full rank. This is 
a very time-consuming task in practice. With the new methods we could obtain this 
information by just looking at the remnant symmetries. We have performed extensive 
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cross-checks in order to show that both methods yield the same results. We have also 
shown how to search for vacua with a given symmetry. Further, we presented the Hibert 
basis method which allows us to survey all D-fiat directions comprised of a selected set 
of fields in very short time. It will be interesting to apply, and to extend, our methods 
to other examples. 
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A Discrete anomalies 

A.l Discrete anomaly calculation 

We calculate the anomaly of the Zf^ symmetry of the configuration discussed in section [31 
The fermions of a superfield with T^f charge q have T^f charge q — 1 because we work in a 
convention where the superpotential carries charge two. Only massless states contribute 
to the anomaly (cf. [21] )• We can hence limit ourselves to the MSSM field content. Since 
all matter fields carry charge one, the corresponding fermion is uncharged and does 
not contribute to the anomaly. Thus, we end up with the following contributions: 



SU(3)c 


SU(2)l 




hu,d 


2-i-3 

C2(3) 


= 3 


gauginos 02(8) = 3 


gauginos 


= 2 


3 


5 



The factor 1/2 in the Higgs contributions is the Dynkin index. The anomaly condition 
is that the sum over all charges be equal mod 2. The total contribution in both cases 
is odd. That is, both symmetries appear anomalous, and the anomaly is, in particular, 
universal, as required for the Green-Schwarz mechanism to work. 

A. 2 Anomaly Mixing 

In our model we have two symmetries which appear anomalous. On the one hand, we 
have the anomalous U(l)anom- On the other hand, the space group selection rule for 
^3 (cf. corresponding to the anomalous space group element in equation (lE.Sp . is 
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anomalous. We wish to answer the question if it is possible to rotate the anomaly 
completely into U(l)anom- 

Consider a setting with U(l) x Zjv symmetry which appear anomalous. We will 
denote the U(l) charge by Q*-*-* and the TLf^ charges by and we will assume all 
charges to be integers. Suppose we have a gauge group G. The anomaly coefficients 
read 

G-G'-U(l) : ^g(^)£(r(^)) = A , (A.la) 
/ 

G-G-Ziv : ^q^^H{r^^^) = B, (A.lb) 
/ 

where the sums run over the irreducible representations of G and i{r) is the Dynkin 
index of the representation r. Specifically, i = 1/2 for the fundamental representation 
of SV{N). 

We can redefine the Tj^ charges by shifting them by integer multiples of the U(l) 
charges. That is, we can define new charges g'^*-* = g'-*-' + nQ^^\ Then the new 
G — G — anomaly coefficient is given by 

= ^ +nQ(^)) £(r(^)) = B + nA (A.2) 
/ / 
with n E 1j. Anomaly freedom requires 



=0 modr^ where V = { ^^^^^ ■ (A.3) 

Hence, the can be made anomaly-free if there is a solution to 

B + nA = mod 1] (A.4) 
for n G Z. 

The anomaly coefficient for the anomalous U(l)anom in our model is given by 

SU(3)c - SU(3)c - U(l)anom : A = 15 (A.5) 

for the SM gauge group SU(3)c. The U(l)anom charges are normalized to integers. The 
anomaly coefficients are the same for SU(2)l and for the gauge group factors of the 
hidden sector because of the Green-Schwarz mechanism. In addition, the anomaly 
coefficient turns out to be 

SU(3)c - SU(3)c - : B = ^. (A.6) 

The equation 

B + nA = ^ + 15?2 = modi (A.7) 

has no solution for n G Z and therefore the anomaly of cannot be removed. Alto- 
gether we have demonstrated that there are (at least) two independent anomalies, i.e. 
the imaginary part of the dilaton shifts both under U(l)anom and 7^2^ transformations 
and there is no way of rotating the Tj^^ anomaly into U(l)anom- 
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B Identifying the symmetry 

In section [3] we explained our strategy for constructing promising vacuum configurations. 
A crucial point is the identification of vacua that exhibit a TL^ symmetry, with charges 
given in table 13.11 In this appendix, we will give some details of how this identification 
is done. 

In a particular vacuum, our models exhibit Abelian discrete symmetries which are 
calculated by the methods described in [30]. A finite Abelian group can always be written 
as a direct product G = Hp^ x . . . x Hp^ where the pi are pairwise distinct primes and 
Hp = Tjpei X ... X Tjpem wlth ci < " " " < Cm positive integers. In a first step, we have to 
make sure that G has a Z4 subgroup. This can be done by looking at the subgroup H2 
of G. 

In the next step, we want to know whether the is of R or non-i? type. In order to 
answer this question, we only have to look at the transformation of the superpotential 
under the Z4; thus one can unambiguously see if the is R or non-i?. 

If the considered vacuum exhibits a symmetry, we have to check that the charges 
of the matter fields match the ones in table 13. 1[ A technical problem appears if H2 
consists of more than one factor, i.e. if there is an additional symmetry where n is 
a positive integer. In this case there are many equivalent charge assignments some of 
which will make the obvious while others will conceal its existence. This freedom 
corresponds to the automorphisms of G [49j whose number can be large, e.g. x Z4 x Z4 
has 1536 automorphisms. An important fact is that the automorphism group factorizes 
in the way we have written G, i.e. Aut(ifpj x Hp^) = Aut(ifpJ x Ant^Hp^). Thus, we 
only need to look at H2. 

The automorphisms can be represented by certain matrices acting on charge vectors 
[19] . To see whether the is present in a vacuum, we scan over all possible charge 
assignments of H2 and look for a Z4 subgroup under which all SM matter fields have 
charge 1. The other states can be even or odd, as long as they are vector-like and can 
be decoupled. 

To illustrate this procedure, let us look at a simple example in table IB. II The two 



charge assignments are equivalent. While in (a) it is not obvious that there is a Z4 
subgroup under which all ip have charge 1, in the second charge assignment in (b) all 
fields have charge 1 with respect to the second Z4 subgroup. 



C Hilbert bases and D-flatness 
C.l General discussion 

In this appendix, a simple method is described that allows us to analyze D-fiatness. 

It is well known that D-fiat directions correspond to holomorphic gauge invariant 
monomials [50l|5T]. We will present now a method to compute all holomorphic gauge 
invariant monomials. Let us look at a theory with gauge group U(l)". As can be 
easily seen, a monomial 0"^02^ . . . 0^'" is D-fiat under the j^^ U(l) factor if (cf. e.g. [161 
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(a) 





Z2 


Z4 


Z4 


^1 


1 





3 




1 


3 


2 


^3 








3 



(b) 





Z2 


Z4 


Z4 







1 


2 




1 


1 


3 


^3 


1 


1 






Table B.l: An example for a hidden Z4 symmetry under which all fields have charge 1. 
The two charge assignments are equivalent. 

appendix B, equation (B3)]) 

J2q^^n. = (C.l) 

i 

where q^''"' is the charge of the field (pi under the j^^ U(l). The index i runs from 1 to the 
number of fields k in the monomial and the index j over the number of different U(l) 
factors. We can rewrite this as a matrix equation, 



(C.2) 



where rij counts how often a field occurs in the monomial. The charges can always 
be scaled to become integers whereas the field multiplicity rii has to be a non-negative 
integer. Thus, the condition to have a D-fiat monomial is 

Q-x = 0, Qe Znxk , X e (C.3) 

This is a system of homogeneous linear Diophantine equation over non-negative integers. 
Such equations can be solved completely by constructing the corresponding Hilbert basis 
(see for example [52]). A Hilbert basis 'H(Q) is a complete set of all minimal solutions 
to equation (lC.3p . A solution is called minimal if it is non-trivial and there exists no 
smaller solution. Here 'smaller' means that, given a solution x, there is no other solution 
y ^ X with Hi < Xi for all i = 1 . . .k. Given ^{Q) we can construct all non-negative 
solutions to equation flC.3|) by forming linear combinations of the basis solutions with 
non- negative integer coefficients. The Hilbert basis ^HiQ) is therefore a basis for all 
D-fiat monomials. In practice the Hilbert basis for a given matrix Q can be computed 
with the help of computer algebra packages like [53]. Such packages can be applied even 
for rather large matrices. 

Let us look at an example from [25] with four fields. The charges of the fields are 
summarized in the charge matrix 

Q = il \ V) ■ 
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The Hilbert basis T-LiQ) is given by the three vectors 



X 



/1\ 





\V 



2 
2 



1 
1 

V/ 



x,y,z G niQ) 



(C.5) 



4^' can be characterized by four- 



All holomorphic gauge invariant monomials 

vectors w = (ni, n2, ns, n^)'^, which are given by w = ax + f3 y + •y z with a, /3, 7 G N. 
We recognize an important property of the Hilbert basis: while the dimension of D- 
flat directions is 2, i.e. the number of fields minus the number of independent D-term 
constraints, the length of the Hilbert basis is larger, namely 3. There is one relation 
between the Hilbert basis elements, 



X + y 



2z 



(C.6) 



This is the price one has to pay for being able to express any D-flat direction as an 
integer linear combination of basis monomials. 



C.2 Hilbert basis for the vacuum discussed in section [3] 

With the Hilbert basis method we could identify a complete set of D-flat directions 
composed of fields in equation fl3.2p . As discussed, the key feature of this Hilbert basis 
is that any gauge invariant holomorphic monomial can be expressed as product of the 
basic monomials. We obtain 6184 monomials. An example of a monomial, which has 
negative charge under the anomalous U(l), is given in equation (ID.9p . 



D Details of the supersymmetric vacuum configura- 
tion 

D.l D-flatness 

Hidden SU(2) breaking. We first look at SU(2). We find that we can switch on 
1/3, . . . , ?/6- This leads to the 6 SU(2) invariant monomials 

{^(*)}SU{2) = {2/32/4,2/32/5,2/32/6,2/42/5,2/42/6,2/52/6}- (D-1) 

There is one relation between the monomials, 

(2/32/4) (2/52/6) - (2/32/5) (2/42/6) + (2/32/6) (2/42/5) = , (D.2) 

such that the number of fiat directions is 5, which is consistent with 8 components being 
switched on whereby 3 directions get eaten by the SU(2) gauge multiplet. The monomial 
yiy2, which also has R charge 0, vanishes. A possible way to have the above monomials 
non- vanishing is to set 

y? = y? = y? = 2/P = y? = y? = o (D.3a) 
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and to have, correspondingly, 

yt\y':},yr\yr' ^ o. (D.sb) 

Hidden SU(3) breaking. There are 16 SU(3) invariant composites of hidden SU(3) 
3- and 3-plets with R charge 0, and can, from this perspective, acquire a VEV, namely 

{^o'-'}sU(3) = {xiXi, X2X3, 0:2X4, X2X5, X3X4, XsXs, X4X4, 0:4X5, 0:5X4, X5X5, 

X1X2X3, X1X3X4, X1X3X5, X2X3X4, X2X3X5, X2X4X5} . (D.4) 

There are different branches of SU(3) flat directions. In what follows we discuss one 
particular of them. 

In this branch there are 13 composites with non-zero VEV, 

{.y#*^'^}su(3) = {xiXi, X2X3, X2X4, X2X5, X3X4, X3X5, X4X4, X4X5, X5X4, X5X5, 

XiX2X3,XiX3X4,XiX3X5} , (D.5) 

while the other 3 monomials vanish, 

(X2X4X5) = (X2X3X5) = (X2X3X4) = . (D.6a) 
Assuming that the above {^^*^}su(3) VEVs do not vanish, we arrive at the relations 

/_ X (X3X5) (X1X3X4) I _ . (X3X5) (X1X2X3) 

(a;3a;4) = 7 ^ , (X2X3) = 7 r , (D.6b) 

_ . (X1X2X3) (X4X5) _ (X1X3X5) (X4X5) 

(X2X5) = r , (X5X5) = r , (D.6c) 

(X1X3X4J (,XiX3X4j 

_ . (X1X2X3) (X4X4) _ (X1X3X5) (X4X4) 

(X2X4) = r , (X5X4) = r . (D.6d) 

(X1X3X4J (X1X3X4J 

This leaves us with 7 independent SU(3) monomials, a possible choice is given by 

{0:3X5, XiXi, X4X4, X4X5, X1X2X3, X1X3X4, X1X3X5} . (D.7) 

So we see explicitly that this branch of D-fiat directions has dimension 7. This is in 
agreement with the result obtained with the STRING VACUA package |34j . 

We can satisfy the constraints by setting (in an appropriate gauge) various compo- 
nents to zero. The only non-vanishing components are 

^[3] JA JA JA J3] [2] [2] [2] , . 

Abelian singlets. We can switch on the Abelian singlets 

{^o'-*} singlet = {01, 02, 03, 04, 05, 06, 07, 08, 09, 010, 011, 012, 013, 014} • 

Since there are 8 U(l) factors that get broken, there are 8 additional D-term constraints. 
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Cancellation of the FI term. With the Hilbert basis method we were able to com- 
pute all gauge invariant monomials carrying negative charge with repect to U(l)anom- 
An example is 

^FI = 0^104 0? 08 09- (D-9) 



D.2 F-flatness 

Remnant symmetry. Switching on the above fields breaks the gauge, R and 
other discrete symmetries down to Gsm x ^4 x [SU(2)]. We decompose the moduli 
space in SU(3) and SU(2) composites and basic fields ^r™"* where r denotes the 
charge of the corresponding objects. A prominent role will be played by the singlet 
fields with R charge 2, which are given by 

{J^t]s\]{i) = {a;iX3,a;iX4,XiX5,XiX2,Xia;4,XiX5,X2X3, 

X2X3X4, X2X3X5, X3X4X5, X1X2X3, X1X2X4, X1X2X5} , 

{^t}sV{2) = {1/12/3, 2/1I/4, Z/12/5, Z/22/3, Z/22/4, 2/2I/5, Z/22/6} , 

{-^2*''}smglet = {01, 02, 03, 04, 05, 06, 07>08> 09, 010, 011, 012} • (D.IO) 



i^'-term constraints. One can use the above monomials for counting the independent 
F-term constraints. As discussed in section [21 the superpotential will be of the form 



(1) 



+ . 



(D.ii; 



where the omission contains only terms which are at least quadratic in the 

ft' are some functions of the monomials with R charge 0. The potentially non-trivial 
F-terms are then 



2-^ rl^. 



ft' 



(D.12) 



2 t \^n-2 /singlet 



as we look at vacua with unbroken Z^, i.e. — 0- For 

each equation gives a non-trivial constraint on the ^q"'^ The number of independent 
constraints is given by the rank of the matrix 



900') 



(D.13) 



where the (p^^'' comprise all component fields appearing in monomials, evaluated at the 
vacuum. 

We evaluated the rank of the matrix for the SU(3) and SU(2) monomials in the 
vacuum defined by (]D.3|) and (ID.Sp . The result is that there are 7 independent F-term 
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constraints in the SU(3) case and 4 in the SU(2) case. Adding the constraints from 
the non-Abehan singlets with R charge 2, we therefore obtain 7 + 4 + 12 = 23 F-term 
conditions. At this point, the supersymmetry conditions seem to over-constrain the 
system, as the number of D-flat directions is 7 + 5 + 14 — 8 = 18. Note, however, that 
there are 6 additional degrees of freedom which we have not discussed yet: the Tj- and 
f/j-moduli of our x orbifold. The functions /g"^'' will also depend on these fields, 
which obviously have R charge (cf. the discussion in [15]). Using these additional 
degrees of freedom we will generically be able to satisfy the constraints, and generically 
there will be only one flat direction in the moduli space formed out of the standard 
model singlet degrees of freedom^ 

There are also SU(3) invariant composites with odd charge. The superpotential 
will contain terms of the form 

W D ^^i(™^-^i(")-/^'"^(^o^'\...) (D.14) 

m,n 

and analogous terms for the This will then lead to non-trivial mass terms for 

the vanishing SU(3) triplets and anti-triplets. Analogous statements hold for the other 
fields with odd R charges. 

In summary, we expect that the vacuum discussed here is such that supersymmetry 
conditions can be satisfied. Moreover, we find that all but one of the fields are fixed by the 
D- and F-term constraints. Unlike in the case without a residual symmetry, due to 
the the superpotential expectation value is guaranteed to vanish at the perturbative 
level. 



E Details of the model 

The orbifold model is defined by a torus lattice that is spanned by six orthogonal vectors 
Co, a = 1, . . . , 6, the x twist vectors Vi = (0, 1/2, -1/2) and V2 = (-1/2, 0, 1/2), 
and the associated shifts 

Vi = (^-^, 0, 0, 0, 0, 0, 0^ (0, 0, 0, 0, 0, 0, 0, 0) , (E.la) 

1 1 
2'~2' 



V2 = (0,^,-^,0,0,0,0,0) (0,0,0,0,0,0,0,0) , (E.lb) 



''Of course, there is the hidden SU(2) sector which contains further massless degrees of freedom. We 
kept this SU(2) unbroken on purpose as it may serve as a toy hidden sector for dynamical supersymmetry 
breakdown. 
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and the six discrete Wilson lines 



Wi = (0^) (0«) , (E.2a) 
/3 111111 1\/ 1111 \ 



4' 4' 4' 4' 4' 4' 4' 4J 4'4' 4' 4' 4' 4' A' Aj ' (E-2c) 
/3 1 3 1 1 3 1 1\ / 3 1 1 3 3 1 3 3\ 

^6 = U.o>-:T'-o.-o'-o.o>o -o.-o>-:T>:T'-o>-:T>-o'o (E-2d) 



^2'2' 2' 2' 2' 2'2'2yV 2' 2' 2'2' 2' 2' 2'2^ 
W2 = = We, (E.2e) 

corresponding to the six torus directions Cq,. Additionally, we divide out the symmetry 
corresponding to 

r = ^(62 + 64 + 66) (E.3) 
with a gauge embedding denoted by W (the freely acting Wilson line) where 

W = ^(1^2 + + 1^6) = ^W2. (E.4) 

The anomalous space group element reads 

fl'anom = {k , £] Ui, n2, u^, n4^, u^, uq) = (0, 0; 0, 0, 1, 0, 0, 0) , (E.5) 
where the boundary conditions of twisted string are 

X(r,cr + 27r) = J X{t, a) + n^e^ (E.6) 

with and u denoting the rotations corresponding to Vi and V2- The spectrum is given 
in table EH In addtion there are 37 GgM x [SU(3) x SU(2) x SU(2)]hid singlets. In 
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Qi 


Ui 
















Vi 
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3 
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3 
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3 


5 
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6 


SU(3)c. 


3 


3 
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1 


1 


SU(2)l 


2 
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1 
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2 
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1 


U(l)y 


1 

6 


2 
3 


1 

3 


1 

3 


1 

2 


1 

2 


1 














SU(3) 


1 


1 


1 


1 


1 


1 


1 


3 


3 


1 


1 


SU(2) 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


1 


SU(2) 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 



Table E.l: The states with their quantum number w.r.t. the SM and the hidden sector. 

table IE. 21 we list the full spectrum. In addition to the states shown there, the spectrum 
contains the following (untwisted) moduli: the dilaton S, three Kahler moduli Tj and 
three complex structure moduli Ui. 
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Table E.2: Spectrum of the model at the orbifold point. The last two columns 
list the (g)eneral and the (c)onfiguration labels. If there are two labels in one 
line, this corresponds to the twist parameter m = for the first label and 

ni = 1 for the second. The two states form a doublet under a D4 symmetry. 
In this model the three charges (corresponding to the three orbifold 
planes) of the respective sectors read: R{Ui) = (2,0,0), R(U2) = (0,2,0), 
R{Ui) = (0,0,2), iJ(r(i,o)) = (0,1,1), R(T(o,i)) = (1,0,1) and ii(T(i,i)) = 
(1,1,0). 
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